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1 The strong CP problem 

The axion was postulated two decades ago JTJ to explain why the strong interactions 
conserve P and CP in spite of the fact that the standard model as a whole violates those 



symmetries. Consider the Lagrangian of QCD: 

1 n 

Cqcd = - I C^(^ + ^[^iD ll g j -{m j q+q Bj + h.c.) 
4 i=i 

The last term is a 4-divergence and hence does not contribute in perturbation theory. 
That term does however contribute through non-perturbative effects associated with 
QCD instantons ||. Such effects can make the physics of QCD depend upon the value 
of 9. Using the Adler-Bell-Jackiw anomaly ||, one can show that 6 dependence must be 
present if none of the current quark masses vanishes. (If this 9 dependence were absent, 
QCD would have a Ua(1) symmetry and would predict the mass of the rj' pseudo-scalar 
meson to be less than y/Sm^ « 240 MeV 0, contrary to observation.) One can further 
show that QCD depends upon 9 only through the combination of parameters: 

9 = 9 — arg(mi, m 2 , . . .m n ) (2) 
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If 9 ^ 0, QCD violates P and CP. The absence of P and CP violations in the strong 
interactions therefore places an upper limit upon 9. The best constraint follows from the 
experimental bound on the neutron electric dipole moment which yields: 9 < 1CT 9 . 

The question then arises: why is 9 so small? In the standard model of particle physics, 
the quark masses originate in the electroweak sector of the theory. This sector must violate 
P and CP to produce the correct weak interaction phenomenology There is no reason in 
the standard model to expect the overall phase of the quark mass matrix to exactly match 
the value of 9 from the QCD sector so that 9 < 10~ 9 . In particular, if CP violation is 
introduced in the manner of Kobayashi and Maskawa J7|, the Yukawa couplings that give 
masses to the quarks are arbitrary complex numbers and hence arg det m q and 9 have no 
reason to take on any special value at all. 

The problem why 9 < 1CT 9 is usually referred to as the "strong CP problem". The 
existence of an axion would solve it. There are other solutions. Setting m u = removes the 
^-dependence of QCD and hence the strong CP problem as well. However, m u = may 
cause problems with the successful current algebra relations among pseudo-scalar meson 
masses. I refer the reader to refs.|8|, |J for recent discussions of the issues involved. Another 
type of solution involves the assumption that CP and/or P is spontaneously broken but 
is otherwise a good symmetry. In this case, 9 is calculable and may be arranged to be 
small fll0|] . Finally, let's emphasize that the strong CP problem need not be solved in the 



low energy theory. Indeed, as Ellis and Gaillard JTT| pointed out, if in the standard model 



9 = near the Planck scale, then 9 ^ 10 9 at the QCD scale. 

Peccei and Quinn |L2| proposed to solve the strong CP problem by postulating the 



existence of a global Upq(1) quasi-symmetry. To do its job, Upq(1) must be a symmetry 
of the theory at the classical (i.e., at the Lagrangian) level, it must be broken explicitly by 
those non-perturbative effects that make the physics of QCD depend upon 9, and finally 
it must be spontaneously broken. The axion |13[ is the quasi-Nambu-Goldstone boson 



associated with the spontaneous breakdown of Upq(1). One can show that, if a Upq(1) 
quasi-symmetry is present, then 

9 = 9 - arg{m x . . .m n ) — , (3) 

J a 

where a(x) is the axion field and f a , called the axion decay constant, is of order the vacuum 
expectation value (VEV) which spontaneously breaks Upq(1). It can further be shown 



14|1 that the non-perturbative effects that make QCD depend upon 9 produce an effective 
potential V(9) whose minimum is at 9 = 0. Thus, by postulating an axion, 9 is allowed 
to relax to zero dynamically and the strong CP problem is solved. 

The properties of the axion can be derived using the methods of current algebra [T3L |T5||. 
The axion mass is given in terms of f a by 

m a ~ 0.6 eV . (4) 

la 

All the axion couplings are inversely proportional to f a . For example, the axion coupling 
to two photons is: 

a a(x) -* -> 

C-arr = —#7 7 E ■ B (5) 

7T Ja 
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where E and B are the electric and magnetic fields, a is the fine structure constant, and g 7 
is a mo del- dependent coefficient of order one. g 7 



0.36 in the DFSZ model 16 



-0.97 in the KSVZ model 



whereas 

A priori the value of f a , and hence that of m a , is 



arbitrary. However, searches for the axion in high energy and nuclear physics experiments 
combined with astrophysical constraints, the latter derived by considering the effect of the 
axion upon the lifetimes of red giants and SN1987a, rule out m a >^ 3 10~ 3 eV [||, 18]. In 
addition, as will be discussed in section II, cosmology places a lower limit on m a of order 
10~ 6 eV by requiring that axions do not overdose the universe. 



2 Dark matter axions 

For small masses, axion production in the early universe is dominated by a novel mechanism 
[T9|j . The crucial point is that the non-perturbative QCD effects that produce the effective 



potential V(9) are strongly suppressed at temperatures high compared to Aqcd pOfl - At 
these high temperatures, the axion is massless and all values of (a(x)) are equally likely. 
At T ~ 1 GeV, the potential V turns on and the axion field starts to oscillate about a 
CP conserving minimum of V. These oscillations do not dissipate into other forms of 
energy because, in the relevant mass range, the axion is too weakly coupled for that to 
happen. The oscillations of the axion field may be described as a fluid of axions. The 
typical momentum of the axions in the fluid is the inverse of the correlation length of the 
axion field. Because that correlation length is of order the horizon, we have p a ~ (10~ 6 
sec)" 1 ~ 10~ 9 eV at T ~ 1 GeV. Afterwards, p a decreases in time as R^ 1 where R is 
the cosmological scale factor. Thus the axion fluid is very cold compared to the ambient 
temperature. 

Let me briefly indicate how the present cosmological energy density of this axion fluid 
is estimated. Let ip(x) be the complex scalar field whose VEV v spontaneously breaks 
Upq(1). At extremely high temperatures, the Upq{\) symmetry is restored. It becomes 
spontaneously broken when the temperature drops below a critical value Tpq of order v. 
Below Tpq, the axion field a(x) appears as the phase of the VEV of ip: (y?(x)) = ve ia ^' v . 
We must now distinguish two cases. Either inflation occurs with reheat temperature 
below Tpq, or not (i.e., inflation does not occur or it occurs with reheat temperature 
above Tpq). In the first case, inflation homogenizes the axion field and there is only 
one contribution to the axion cosmological energy density, the contribution from so-called 
"vacuum misalignment" . In the second case, there are additional contributions from axion 
string and axion domain wall decay. Only the contribution from vacuum misalignment is 
discussed in any detail here. 

When the axion mass turns on near the QCD phase transition, the axion field starts to 
oscillate about one of the CP conserving minima of the effective potential. The oscillation 
begins approximately at cosmological time t 1 such that tim a (T(ti)) = 0(1), where m a {T) is 
the temperature dependent axion mass. Soon after t\, the axion mass changes sufficiently 
slowly that the total number of axions in the oscillations of the axion field is an adiabatic 
invariant. T\ = T(t\) has been estimated to be of order 1 GeV. The number density of 
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axions at time t\ is 

n a {h) ~ ^maih^a 2 ^)) ~ 7r/ a 2 i (6) 

where / a = is the axion decay constant introduced earlier. iV is an integer which 
expresses the color anomaly of Upq(1). N also equals the number of CP conserving 
vacua [pi] at the bottom of the 'Mexican hat' potential, i.e., in the interval < - < 2tt. 
In Eq. (6), we used the fact that the axion field a(x) is approximately homogeneous on the 
horizon scale t±. Wiggles in a(x) which entered the horizon long before t\ have been red- 



shifted away [ 22[ . We also used the fact that the initial departure of a(x) from the nearest 
minimum is of order ~ = f a . The axions of Eq. (6) are decoupled and non-relativistic. 
Assuming that the ratio of the axion number density to the entropy density is constant 



from time ti till today, one finds |19 



/0.6 10" 5 eV\^ /200 MeV\ * f 75 km/s-Mpc\ 2 

a - = {—^—) l^rj { — k — J (7) 

for the ratio of the axion energy density to the critical density for closing the universe. Hq 
is the present Hubble rate. Eq. (7) implies the bound m a >^ 1CT 6 eV. 

However, there are many sources of uncertainty in the estimate of Q a . The axion energy 
density may be diluted by the entropy release from heavy particles which decouple before 
the QCD epoch but decay afterwards P5[ , or by the entropy release associated with a first 
order QCD phase transition. On the other hand, if the QCD phase transition is first order 



24]] , an abrupt change of the axion mass at the transition may increase f2 a . If inflation 
occurs with reheat temperature less than T P q, there may be an accidental suppression of 
Q a because the homogenized axion field happens to lie close to a CP conserving minimum. 
Because the RHS of Eq. (7) is multiplied in this case by a factor of order the square of the 
initial vacuum misalignment angle which is randomly chosen between — tt and +n, 

the probability that Q a is suppressed by a factor x is of order yfx. This rule cannot be 
extended to arbitrarily small x however because quantum mechanical fluctuations in the 



axion field during the epoch of inflation do not allow the suppression to be perfect |25 



If inflation occurs with reheating temperature larger than Tpq or if there is no inflation, 



there are contributions to fl a from axion string 26 and axion domain wall decay E7[ in 



addition to the contribution, Eq. (7), from vacuum misalignment. My collaborators and 



I [PHI have estimated each of these additional contributions to be of the same order of 



magnitude as that from vacuum misalignment. Other authors [26|, 2|| have estimated that 



the contribution from axion string decay dominates over that from vacuum misalignment 
by a factor 100 or a factor 10. 

The axions produced when the axion mass turns on during the QCD phase transition 
are cold dark matter (CDM) because the axions are non-relativistic from the moment of 
their first appearance at 1 GeV temperature. Studies of large scale structure formation 
support the view that the dominant fraction of dark matter is CDM [[JI|. Moreover any 



form of CDM necessarily contributes to galactic halos by falling into the gravitational wells 
of galaxies. Hence, there is excellent motivation to look for CDM candidates as constituent 
particles of our galactic halo, even after some fraction of our halo has been demonstrated 
to be in MACHOs flHTIl or some other form. 
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Finally, let's mention a particular kind of dumpiness 132 1 which affects axion dark 
matter if there is no inflation after the Peccei-Quinn phase transition. This is due to 
the fact that the dark matter axions are inhomogeneous with 5p/p ~ 1 over the horizon 
scale at temperature T\ ~ 1 GeV, when they are produced at the start of the QCD 
phase-transition, combined with the fact that their velocities are so small that they do 
not erase these inhomogeneities by free-streaming before the time t eq of equality between 
the matter and radiation energy densities when matter perturbations can start to grow. 
These particular inhomogeneities in the axion dark matter are immediately in the non- 
linear regime after time t eq and thus form clumps, called 'axion mini-clusters' J32|. These 





have mass M mr ~ 10 13 Mm and size l mr — 10 12 cm. 



3 The cavity detector of galactic halo axions 

Axions can be detected by stimulating their conversion to photons in a strong magnetic 
field [p3| . The relevant coupling is given in Eq. (5). In particular, an electromagnetic cavity 
permeated by a strong static magnetic field can be used to detect galactic halo axions. The 
latter have velocities (3 of order 10~ 3 and hence their energies E a = m a + |m a /3 2 have a 
spread of order 10~ 6 above the axion mass. When the frequency uj = 2k f of a cavity mode 
equals m a , galactic halo axions convert resonantly into quanta of excitation (photons) of 
that cavity mode. The power from axion — > photon conversion on resonance is found to 
be 132 



P = {- 9 -A VB 2 oPa C—Mm(Q L ,Q a ) 
0.5 l(T 26 Watt ( V ^ B ° 



,500 liter/ V7 Tesla 



where V is the volume of the cavity, B is the magnetic field strength, Q L is its loaded 
quality factor, Q a = 10 6 is the 'quality factor' of the galactic halo axion signal (i.e. the 
ratio of their energy to their energy spread), p a is the density of galactic halo axions on 
Earth, and C is a mode dependent form factor given by 



C 



., - -2 



B 2 VJ v d3xe\E w \2 



(9) 



where B (x) is the static magnetic field, E w {x)e lult is the oscillating electric field and e is 
the dielectric constant. 

Because the axion mass is only known in order of magnitude at best, the cavity must 
be tunable and a large range of frequencies must be explored seeking a signal. The cavity 
can be tuned by moving a dielectric rod or metal post inside it. Using Eq. (8), one finds 
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Figure 1: Regions in mass - (coupling/mass) 2 space which have been ruled out by the RBF 
and UF experiments (hatched), which the LLNL experiment has searched so far (densely 
hatched) and which it expects to rule out in the future (shaded border). g aTy = The 
DFSZ and KSVZ model predictions are also shown. 



that to perform a search with signal to noise ratio s/n, the scanning rate is: 
df 12GHz [An 




(10) 

where T n is the sum of the physical temperature of the cavity plus the noise temperature of 
the microwave receiver that detects the photons from a — > 7 conversion. Eq. (10) assumes 
that Ql < Q a and that some strategies have been followed which optimize the search rate. 
The best quality factors attainable at present, using oxygen free copper, are of order 10 5 
in the GHz range. To make the cavity of superconducting material is probably not useful 
since it is permeated by a strong magnetic field in the experiment. 

Eq. (8) shows that a galactic halo search at the required sensitivity is feasible with 
presently available technology, provided the form factor C can be kept at values of order one 
for a wide range of frequencies. For a cylindrical cavity and a longitudinal magnetic field, 
C = 0.69 for the lowest TM mode. The form factors of the other modes are much smaller. 
The resonant frequency of the lowest TM mode of a cylindrical cavity is /= 115 MHz (jjf) 
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where R is the radius of the cavity. Since 10" 6 eV = 2tc (242 MHz), a large cylindrical 
cavity is convenient for searching the low frequency end of the range of interest. To extend 
the search to high frequencies without sacrifice in volume, one may power-combine many 
identical cavities which fill up the available volume inside a magnet's bore |35| , |36| . This 
method allows one to maintain C = 0(1) at high frequencies, albeit at the cost of increasing 
engineering complexity as the frequency, and hence the number of cavities, is increased. 

Pilot experiments were carried out at Brookhaven National Laboratory |37] and at 
the University of Florida |38[ which demonstrated the feasibility of the cavity detection 
method. The (magnetic field) 2 x volume provided by the magnets used in these experi- 
ments were relatively low: B^V = 0.36T 2 m 3 and 0.45T 2 m 3 respectively for Brookhaven 
(RBF) and Florida (UF). Fig. 1 shows the limits that these experiments placed on the 
coupling g ail = ^j- as a function of the axion mass m a assuming that the local density 

of galactic halo axions is p a = | ■ i0~ 24 gr/cm 3 . 

Second generation experiments are presently under way at Lawrence Livermore Na- 
tional Laboratory (LLNL)[]39| and at Kyoto University [[5J. The LLNL experiment uses a 
much larger magnet {BqV = 12T 2 m 3 ) than the pilot experiments. It started taking data 
in Feb. '96. Up till now, it has searched the range of frequencies 670 to 800 MHz, or 2.77 
to 3.30 /ieV (see Fig. 1), at a level of sensitivity sufficient to discover KSVZ axions if they 
are the constituents of our halo. The LLNL experiment plans to cover the mass range 
1.3 < m a < 13yueV. For the high frequency end of this range, it will use an array of four 
cavities kept in tune and whose outputs are power-combined. SQUIDs are being devel- 
opped by members of the collaboration as front-end detectors of the microwave power from 
axion to photon conversion. So far, the LLNL experiment has used the same microwave 
detection technology, HEMT amplifiers, as the pilot experiments. SQUIDs may boost the 
sensitivity by a factor of ten or more. This would allow the LLNL experiment to search 
for DFSZ axions whereas only sensitivity to KSVZ axions had been originally planned. 

The Kyoto experiment has a magnet of size similar to that of the pilot experiments 
but uses a beam of Rydberg atoms to count the photons from a — > 7 conversion. Single 
photon counting constitutes a dramatic improvement in microwave detection sensitivity. 
With HEMT amplifiers one needs to have thousands of a — > 7 conversions per second and 
integrate for about 100 sec to find a signal in the noise. With single photon counting, a 
few a — > 7 conversions suffice in principle. To build a beam of Rydberg atoms capable 
of single photon counting is a considerable achievement in itself. In addition, the cavity 
will be cooled by a dilution refrigerator down to a temperature (~ 10 mK) where the 
thermal photon background is negligible. The Kyoto experiment will first search near 
m a = 10~ 5 eV. Its projected sensitivity is sufficient to discover DFSZ axions even if their 
local density is only | of the local halo density. 



4 Caustic rings in the density distribution of cold 
dark matter halos 

If a signal is found in the cavity detector of galactic halo axions, it will be possible to 
measure their energy spectrum with great precision and resolution. Hence there is good 



7 



motivation to ask what can be learned about our galaxy from analyzing the signal. Con- 
versely, if the energy spectrum of cold dark matter particles on Earth can be deduced 
from studies of galactic halo formation, that information may be used to increase the 
signal/noise ratio of the search. It may also be useful in searching for WIMPs. Indeed, 
because they have moved purely under the influence of gravity since their decoupling, 
all cold dark matter particles have exactly the same phase space distribution and hence 
the same energy spectrum on Earth in the limit where their, at any rate tiny, primordial 
velocity dispersions are neglected. 

In many past discussions of dark matter detection on Earth, it has been assumed that 
the dark matter particles have an isothermal distribution. Thermalization has been argued 
to be the result of a period of "violent relaxation" following the collapse of the protogalaxy. 
If it is strictly true that the velocity distribution of dark matter particles is isothermal, 
which seems to be a strong assumption, then the only information that can be gained 
from its observation is the corresponding virial velocity and our own velocity relative to 
its standard of rest. If, on the other hand, thermalization is incomplete, a signal in a dark 
matter detector may yield additional information. 



J.R. lpser and I discussed EJJ the extent to which the phase-space distribution of 
cold dark matter particles is thermalized in a galactic halo and concluded that there are 
substantial deviations from a thermal distribution in that the highest energy particles have 
discrete values of velocity. There is one velocity peak on Earth due to dark matter particles 
falling onto the galaxy for the first time, one peak due to particles falling out of the galaxy 
for the first time, one peak due to particles falling into the galaxy for the second time, etc. 
The peaks due to particles that have fallen in and out of the galaxy a large number of 
times in the past are washed out because of scattering in the gravitational wells of stars, 
globular clusters and large molecular clouds. But the peaks due to particles which have 
fallen in and out of the galaxy only a small number of times in the past are not washed 
out. 

I. Tkachev, Y. Wang and I have used the self-similar infall model of galactic halo 
formation to estimate the local densities and the velocity magnitudes of the dark matter 



particles in the velocity peaks El. We generalized the existing version of that model to 



take account of the angular momentum of the dark matter particles. In the absence of 
angular momentum, the model produces flat rotation curves for a large range of values of 
a parameter e. We find that the presence of angular momentum produces an effective core 
radius, i.e., it makes the contribution of the halo to the rotation curve go to zero at zero 
radius. The presence of a core radius is consistent with observation. The model provides 
a detailed description of the large scale properties of galactic halos including their density 
profiles, their extent and their total mass. By fitting the model to our galactic halo, we 
obtained the average values of the local densities and the velocity magnitudes of the dark 
matter particles in the velocity peaks on Earth. The averages are over all locations at the 
same distance (8.5 kpc) from the galactic center as we are. 

In the course of this study, it was found that caustic rings EJ, Q form in the coherent 
flows of dark matter particles associated with the velocity peaks. A caustic is a location in 
physical space where the density is enhanced because the 3-dim. sheet on which the dark 
matter particles lie in 6-dim. phase space folds back there. The Zel'dovich "pancake" is 
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an example of caustic. It is a surface in physical space where the density diverges as 
x being the distance to the surface. The ring caustic is a closed line, or loop, where the 
density diverges as ~, x being the distance to the line. The line may however have internal 
structure. This remains to be investigated. 

A caustic ring is associated with each flow of dark matter particles in and out of the 
central parts of the galaxy. Thus, there is a caustic ring due to particles falling through 
the galaxy for the first time, a caustic ring of smaller size due to particles falling through 
for the second time, a yet smaller ring due to particles falling through for the third time, 
and so on. For an arbitrary angular momentum distribution of the infalling particles, the 
ring is a closed line of arbitrary shape. However, if the angular momentum distribution 
is dominated by a smooth component which carries net angular momentum, the ring 
resembles a circle. If there is no angular momentum at all, the ring reduces to a point 
at the galactic center. The caustic rings get smeared by the velocity dispersion of the 
infalling dark matter. However, for a galaxy like our own, that velocity dispersion has to 
reach a few 10 km/s for the caustic rings to be completely washed out. 

The self-similar infall with angular momentum model of galactic halo formation predicts 
the radii of the caustic rings in terms of the rotation velocity v rot of the galaxy and the 
model parameters e, j max and h. In theories of large scale structure formation with cold 
dark matter and a Harrison-Zel'dovich spectrum of primordial density perturbations, e is 
in the range 0.2 to 0.3. j max is the maximum value of the dimensionless angular momentum 
distribution of the infalling particles, as defined in ref. f£J. h is the Hubble constant 
in units of 100 km/s.Mpc. It turns out that, in the range e ~ 0.2 to 0.3, the ratios a n /a n _i 
of successive ring radii are nearly independent of v rat , e, h and j ma x- They have the values 

{a„/a n _! : n = 2, 3, 4, 5...} = (0.49, 0.67, 0.76, 0.81, ...) (11) 

which are therefore predictions of the self-similar infall with angular momentum model. 
It is natural to expect the caustic rings of a spiral galaxy to lie in its galactic plane. In 
that case, the caustic rings cause bumps in the galactic rotation curve which are large 
enough to be observed. It turns out that the rotation curve of NGC3198, one of the best 
measured and often cited as providing compelling evidence for halos of dark matter, has 
three bumps at radii whose ratios match the first two entries on the RHS of Eq. (11). 
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